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Abstract 

We study the effect on leptogenesis due to B — L cosmic strings of a 17(1) b~l extension of the 
Standard Model. The disappearance of closed loops of B — L cosmic strings can produce heavy right 
handed neutrinos, N} j’s, whose CP-asymmetric decay in out-of-thermal equilibrium condition can 
give rise to a net lepton (L) asymmetry which is then converted, due to sphaleron transitions, to a 
Baryon ( B ) asymmetry. This is studied by using the relevant Boltzmann equations and including 
the effects of both thermal and string generated non-thermal Nr s. We explore the parameter 
region spanned by the effective light neutrino mass parameter rh\, the mass M\ of the lightest 
of the heavy right-handed neutrinos (or equivalently the Yukawa coupling h±) and the scale of 
B — L symmetry breaking, tjb-l, and show that there exist ranges of values of these parameters, 
in particular with gs-L > 10 11 GeV and h\ > 0.01, for which the cosmic string generated non- 
thermal N} j’s can give the dominant contribution to, and indeed produce, the observed Baryon 
Asymmetry of the Universe when the purely thermal leptogenesis mechanism is not sufficient. We 
also discuss how, depending on the values of gs-L, m\ and h\, our results lead to upper bounds 
on sin 5, where 5 is the the CP violating phase that determines the CP asymmetry in the decay 
of the heavy right handed neutrino responsible for generating the L-asymmetry. 
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I. INTRODUCTION 


], A 0 


Present low energy neutrino oscillation data in, a a are elegantly explained by the 
neutrino oscillation hypothesis with very small masses (< 1 eV) of the light neutrinos. 
Neutrinos can have either Dirac or Majorana masses. Small Majorana masses of the ligh 


neutrinos, however, can be generated in a natural way through the seesaw mechanism [4] 
without any fine tuning. This can be achieved by introducing right handed neutrinos, Nrs, 
into the electroweak model which are neutral under the known gauge symmetries. The 
Majorana masses of these Nrs are free parameters of the model and are expected to be 


either at TeV scale |5( or at a higher scale |g, 7]. This indicates the existence of new physics 
beyond Standard Model (SM) at some predictable high energy scale. 

The heavy right handed Majorana neutrinos are also an essential ingredient in currently 
one of the most favored scenarios of origin of the observed barvon asymmetry of the Universe 
(BAU), namely, the “baryogenesis via leptogenesis” scenario j|,P,|lO|, 11]. Majorana mass of 
the neutrino violates lepton number (L) and thus provides a natural mechanism of generating 
a lepton asymmetry in the Universe. Specifically, leptogenesis can occur via the L- violating, 


CP-asymmetric, out-of-equilibrium j3] decay of the Nr s into SM leptons and Higgs. The 
resulting L-asymmetry is then partially converted to a baryon (P)-asynnnetry via the non 
perturbative B + L violating (but B — L conserving) electroweak sphaleron transitions jl^ . 
The attractive aspect of the leptogenesis mechanism is the link it implies between the physics 
of the heavy right handed neutrino sector and the experimental data on light neutrino flavor 
oscillation, thus making the scenario subject to experimental tests. Indeed, the magnitude 
of the L (and thus B) asymmetry produced depends on, among other parameters, the 
masses of the heavy neutrinos, which in turn are related to the light neutrino masses via the 
seesaw mechanism. The mass-square differences amongst the three light neutrino species 
inferred from the results of neutrino experiments, therefore, place stringent constraints on 
the leptogenesis hypothesis. 

A natural way to implement the leptogenesis scenario is to extend the SM to a gauge group 
which includes B — L as a gauge charge. The heavy neutrino masses are then determined by 
the scale of spontaneous breaking of this gauge symmetry. Further, with B — L a conserved 
gauge charge and B + L anomalous, we can start with a net B = L = 0 at a sufficiently 
early stage in the Universe. The observed B asymmetry must then be generated only after 





the phase transition breaking the B — L gauge symmetry. 

It is well known that phase transitions associated with spontaneous breaking of gauge 
symmetries in the early Universe can, depending on the structure of the symmetry group and 
its breaking pattern, lead to formation of cosmic topological defects jlj, ll5| of various types. 
In particular, the simplest choice for the B — L gauge symmetry group being a U (1) b-l, the 
phase transition associated with spontaneous breaking of this U(1)b-l hr the early Universe 
would, under very general conditions, lead to formation of cosmic strings jlj, LL5J carrying 
quantized B — L magnetic flux. These U B — L” cosmic strings can be a non-thermal source 
of Nr s whose decay can give an extra contribution to the L and thereby B asymmetry in 
addition to that from the decay of Nrs of purely thermal origin (“thermal” leptogenesis). 
This can happen in the following two ways: 

First, since the Higgs held defining the B—L cosmic string is the same Higgs that also gives 
mass to the N R through Yukawa coupling, the N R s can be trapped inside the B — L cosmic 
strings as fermionic zero modes [3]. Existence of zero energy solutions of fermions coupled 
to a Higgs held that is in a topological vortex string configuration is well-known 12, 12j, and 
has been studied in a variety of models allowing cosmic strings |l£|, [20, l2ll> 122, 123] • As closed 
loops of B — L cosmic strings oscillate, they lose energy due to emission of gravitational 
radiation and shrink in size. Eventually, when the size of the loop becomes of the order of 
the width of the string, the string loop disappears into massive particles among which will 
be the Nrs which were trapped inside the string as zero modes. Each closed loop would be 
expected to release at least one N R when it finally disappears, and the decay, of these Nr s 
would then give a contribution to the BAU through the leptogenesis route |lfil |. 

Second, collapsing, decaying or repeatedly self-intersecting closed loops of cosmic strings 
would in general produce heavy gauge and Higgs bosons of the underlying spontaneously 
broken gauge theory. In the context of cosmic strings in Grand Unified Theories (GUTs), 


H 


iffsrs bosons released from 


25, |2f 


The sphaleron 


the CP asymmetric ^-violating decay of the heavy gauge and 
cosmic string loops would produce a net B asymmetry [ 
transitions would of course erase the B asymmetry unless a net B — L was generated. If 
the strings under consideration are the B — L cosmic strings, which can be formed at an 
intermediate stage of symmetry breaking in a GUT model based on S'O(IO), for example, 
then the heavy gauge and Higgs bosons released from the decaying or collapsing B — L 
cosmic string loops can themselves decay to Nrs since the Nr s have Yukawa and gauge 






couplings to the (string-forming) Higgs and gauge boson, respectively. The decay of these 
Higgs and gauge boson generated Nr' s can produce a net B — L and thus contribute to the 
BAU through the leptogenesis route irrespective of the existence of zero modes of Nr's on 
cosmic strings. 


In a previous work 


28], we made a general analytical estimate of the contribution of the 


non-thermal Nr s produced by B — L cosmic string loops to BAU. It was shown there that, 
in order for the resulting B asymmetry not to exceed the measured BAU inferred from the 
WMAP results j2t|, the mass M\ of lightest right handed Majorana neutrino had to satisfy 


the constraint 


M 1 < 2.4 x 10 12 (?7 B _ L /10 13 GeV) 1/2 GeV, 


( 1 ) 


where t/b-l is the U(1)b-l symmetry breaking scale. In the above mentioned analytical 
study we had (a) taken the CP asymmetry parameter e\ to have its maximum value (see 
below), (b) not taken into account the contribution to BAU from the decay of the already 
existing thermal Nr's, and (c) neglected all wash-out effects (see below) on the final B — L 
asymmetry. Indeed, our aim there was to make a simple analytical estimate of the possible 
maximum contribution of the cosmic string generated non-thermal Nr s to the measured 
BAU. It is, of course, clear that a complete analysis of leptogenesis in presence of cosmic 
strings can only be done by solving the full Boltzmann equations *9, 11] that include the 
non-thermal Nr s of cosmic string origin as well as the already existing thermal Nr' s and 
take into account all the relevant interaction processes including the wash out effects. This 
is the study taken up in this paper. 

The main results obtained in the present paper can be summarized as follows: First, 
we confirm the result, obtained earlier in our analytical study j^, that B — L cosmic 
string loops can give significant contribution to BAU only for t/b-l ~ 10 11 GeV. Second, 
the numerical solution of the relevant Boltzmann equations in the present paper has enabled 
us to track the dynamical evolution of the contribution of the cosmic string generated non- 
thermal Nr' s to the final BAU. Specifically, we find that for sufficiently large values of t/b-l 
and for appropriate ranges of values of the other relevant parameters rhi and h 1 as detailed 
in sec. EVDl the effect of the cosmic string generated non-thermal Nr' s is to produce a 
late-time increase of the final value of BAU as compared to its value in absence of cosmic 
strings. This can be understood from the fact that while the thermal abundance of Nr 's 
decreases exponentially with decreasing temperature, the density of cosmic string generated 





Nrs has a power law dependence (on temperature) inherited from the scaling behavior of 
the evolution of cosmic strings, leading to a domination of the string generated Nr’s over 
the thermal Nr’s at late times for sufficiently large values of t)b-l- hi such situations, we 
are required to place an upper bound on the magnitude of sin 5, where <5 is the CP violating 
phase that determines the CP asymmetry in the heavy neutrino sector, in order not to 
overproduce the BAU. 

The rest of our paper is organized as follows. In section m we briefly review the standard 
thermal baryogenesis via leptogenesis scenario, and discuss the required lower bounds on 
the mass of the lightest right handed neutrino and the B — L symmetry breaking scale. 
In section mu we introduce closed loops of B — L cosmic strings as non-thermal sources of 
Nr s and write down the injection rates of these Nr s due to the two main processes of 
disappearance of cosmic string loops. Section |TV] is devoted to setting up and then solving 
the relevant Boltzmann equations for the evolution of the B — L asymmetry, including the 
non-thermal Nr s of cosmic string origin in addition to the usual thermal ones. The effects 
of the cosmic string generated non-thermal Nr's on the evolution of the B — L asymmetry 
are discussed. Finally, section IVl summarizes our main results. 


II. A BRIEF REVIEW OF BARYOGENESIS VIA LEPTOGENESIS 


A. The general framework 

In the model we consider, the lepton asymmetry arises through the decay of right handed 
Majorana neutrinos to the SM leptons (£) and Higgs (</>) via the Yukawa coupling 

£y = fijU4>N Rj + h.c., (2) 

where f l3 is the Yukawa coupling matrix, and i, j = 1, 2, 3 for three flavors. 

We use the basis in which the mass matrix of the heavy Majorana neutrinos, M, is 
diagonal, where the Majorana neutrinos are defined by Nj = ^=( Nr j + Nr.). However, in 
this basis the Dirac mass matrix mo (= fv, v being the vacuum expectation value of the 
SM Higgs) of the neutrinos is not diagonal. The canonical see-saw mechanism then gives 
the corresponding light neutrino mass matrix 


m v = 


—m£)M 1 m 1 R ) . 


(3) 


In this basis the light neutrino mass matrix m u can be diagonalized by the lepton mixing 
matrix H Ul to give the three light neutrino masses, 


diag(mi,m 2 ,m 3 ) = Ulm v U L - 


( 4 ) 


The best fit values of the mass-square differences and mixing angles of the light neutrinos 
from a global three neutrino flavors oscillation analysis are [31] 


00 = 0,-34°, 6 atm = 0 2 = 45°, 0 3 <13°, 


( 5 ) 


and 


A = \m\ 
A ™aim = l™3 


m , 


7.1 x 1CT 5 eV 2 , 

( 6 ) 

2.6 x 1CT 3 eV 2 . 

(7) 


In the absence of data about the overall scale of neutrino masses, we shall assume throughout 
this paper a normal hierarchy, m, <C m 2 -C m 3 , for the light neutrino masses. With this 
assumption equations © and © give 

m 2 ~ \J Am| = 0.008 eV, 
m 3 ~ v 7 = 0.05eV. 


( 8 ) 


We will use the above value of 7713 to obtain the upper bound on UP-asymmetry in subsection 

©s 


B. Decay of heavy Majorana neutrino and CP-asymmetry 

The C P-asymmetry parameter in the decay of A* is dehned as 

— 


( 9 ) 


r(jv ( et) + r(Ni — (4, t) 

We assume a normal mass hierarchy in the heavy Majorana neutrino sector, M\ < M 2 < M 3 , 
and further assume that the final lepton asymmetry is produced mainly by the decays of 
the lightest right handed Majorana neutrino, A,. The latter is justified |ol llll| because any 
asymmetry produced by the decay of iV 2 and A 3 is erased by the lepton number violating 
interactions mediated by A,. At an epoch of temperature T > M,, all the lepton number 
violating processes mediated by A, are in thermal equilibrium. As the Universe expands and 










cools toL < Mi, the L-violating scatterings mediated by Ah freeze out, thus providing the 
out-of-equilibrium situation |u]] necessary for the survival of any net L-asymmetry generated 
by the decay of the Ah’s. The hnal L-asymmetry is, therefore, given essentially by the 
product of the equilibrium number density of the Ah’s at L ~ Mi and the CP-asymmetry 
parameter e\. The latter is given by HQ 

3 


l e i| — 


-Mimj, sin 8 , 


( 10 ) 


16itv 2 

where v ~ 174 GeV is the electroweak symmetry breaking scale and 8 is the CP violating 
phase. Using the best fit value of m 3 from equation (jHJ) the value of ei can be written as 

Mi \ ( m3 


ei < 9.86 x 1CT 


10 9 GeV / \0.05 eV 


( 11 ) 


C. Thermal leptogenesis and bound on M\ 


The L-asymmetry, created by the decays of Ah, is partially converted to a P-asymmetry 
by the nonperturbative sphaleron transitions which violate B + L but preserve B — L. 
Assuming that sphaleron transitions are ineffective at temperatures below the electroweak 
transition temperature (Lew) 5 the P-asymmetry is related to L-asymmetry by the rela¬ 
tion 


P = 


P 


p — 1 


L ~ -0.55L, 


( 12 ) 


where we have taken p = 28/79 appropriate for the particle content in SM. The net baryon 
asymmetry of the Universe, defined as Y B = ( n B /s ), can thus be written as 


Y b ~ 0.5561^1 T~M 1 d , 

where the factor 0.55 in front comes from equation m, 

v _ n m 

Uvi =-, 


(13) 


(14) 


rijVj being the number density of A^ 1; and d is the dilution factor due to wash-out effects. In 
the above equations s stands for the entropy density and is given by 


2tt 2 

s = — g*T 3 ~ 43.86(^/100)L 3 , 
45 


(15) 


where g* is the total number of relativistic degrees of freedom contributing to entropy of the 
Universe. 










The present-day observed value of the baryon-to-photon ratio rj = (n B — n B )/n 1 inferred 
from the WMAP data is 


WMAP 
'/0 




7.0 V B U 


= (6.118:1) x 10' w . 


(16) 

(17) 


Using equation m in equation m and comparing with (Hzd we get 

1.6 x 1(T 3 ' 


Mi > 10 y GeV x 


Yr 


N i |T~Mi 


0.05eV^ 


m 3 


(18) 


Assuming that N\s are initially (i.e., at T > Mi) in thermal equilibrium and they remain 
so till T ~ Mi, one has Y Nl \ T ~ Ml ~ 2.3 x 1CT 3 x (100/g*). Equation (fTHli then indicates 
that the B — L symmetry breaking scale must satisfy the constraint 7] B _ L >, O(10 9 ) GeV for 
successful thermal leptogenesis, unless we allow the Majorana neutrino Yukawa coupling of 
the lightest right handed neutrino, h\ , to be greater than unity. 


III. B — L COSMIC STRINGS AS SOURCES OF NON-THERMAL HEAVY NEU¬ 
TRINOS 


A. B — L cosmic strings and zero modes of Nr 


The possibility of B — L cosmic strings in various extensions of the SM were studied in 


|34j. For simplicity, we consider in the present paper a model based on the gauge group 
SU(2)i, 0 [/(l)y 0 C/(1)y', where Y' is a linear combination of Y and B — L |35|. We then 
follow the symmetry breaking scheme 


SU(2) l 0 U{1) Y ® U (l)y' (y) = rjB-L SU( 2) L 0 U{ 1) Y 

- > 

(4>) = v U(1 )em, (19) 

- > 

where x i s lh e Higgs boson required to break the U(1)y' gauge symmetry and (j) T = (0 + , 0°) 
is the SM Higgs. 

As the Universe cools below the critical temperature, T B -l , of the B — L symmetry 
breaking phase transition, the Higgs field x develops a vacuum expectation value (x) = Vb-l- 
The same Higgs field also forms strings. The mass per unit length of the string, / x , is given 
by Q fi ~ t)b_ l ~ Tl_ L . The exact value of /i depends on the values of the parameters 







of the model, in particular the relevant gauge and Higgs boson masses, and can differ from 
rj 2 B _ L by up to an order of magnitude depending on the model 36]. Measurements of 
CMB anisotropies have been used to place upper bounds on the fundamental cosmic string 


parameter u in a variety of different cosmic string models; see, for example, 


37]. A recent 


analysis |38| of data on CMB anisotropies and large scale structure together puts the bound 
G/i < 3.4(5) x 10” 7 at 68 (95)% c.l. This most likely rules out cosmic string formation at 
a typical GUT scale ~ 10 16 GeV. However, lighter cosmic strings arising from symmetry 
breaking at lower scales, for example, the B — L cosmic strings discussed in this paper, are 
not ruled out. 

The right handed neutrinos acquire Majorana masses from coupling to the same Higgs 
field x, 

-C x - Nr = i [ihiixN Ri N c Ri + h.c] , (20) 

where h is the Yukawa coupling matrix, and = icr' 2 N B defines the Dirac charge conju¬ 
gation operation. The equation of motion of Nr in the background of y forming the string 
admits \n\ normalizable zero-modes in the winding number sector n ljj, Ujj, Ilia, |22| • On a 
straight string these modes are massless. However, on a generic string they are expected to 
acquire effective masses proportional to the inverse radius of the string curvature. As soon 
as this mass becomes comparable to the mass of the free neutrinos in the bulk medium, 
these neutrinos can be emitted from the string. In any case, when the string loop shrinks 
and finally decays it emits various massive particles: the gauge bosons, the heavy Higgs (y) 
and the massive right handed Majorana neutrinos (A)). 


B. NNs from closed loops of B — L cosmic strings 

A key physical process that governs the evolution of cosmic string networks in an expand¬ 
ing Universe is the formation of sub-horizon sized closed loops which are pinched off from 
the network whenever a string segment curves over into a loop and intersects itself. It is this 
process that allows the string network to reach a scaling regime, in which the energy density 
of the string network scales as a fixed fraction of the radiation or matter energy density in 
the Universe; see, e.g., Ref. jl^ for a text book discussion of evolution of cosmic strings in 
the Universe. 

After their formation, the closed loops eventually disappear through either of the following 






two processes (see Refs. (28, Q for the relevant details): 


1. Slow death 


Closed loops born in non-selfintersecting configurations oscillate freely with oscillation 
time period Lj 2 (L being the length of the loop). In doing so they slowly lose energy due to 
emission of gravitational radiation, and thus shrink in size. Eventually, when the radius of 
the loop becomes of the order of the width (~ of the string — which happens over a 

time scale large compared to L (hence “slow”) — the resulting “tiny” loop loses topological 
stability and decays into elementary particle quanta including the gauge- and Higgs bosons 
associated with the underlying broken symmetry as well the heavy neutrinos iVj’s coupled to 
the gauge- and Higgs bosons. While we expect the final number Nn of the heavy neutrinos 
released per tiny loop to be of order one, it is difficult to be more precise, and we keep this 
number as an undetermined parameter. Assuming that the energy of a loop going into Nrs 
is predominantly in the form of the lightest of the heavy neutrinos, Ni s, the number of 
heavy Majorana neutrinos released from the closed loops disappearing through this “slow 
death” (SD) process per unit time per unit volume at any time t (in the radiation dominated 
epoch) can be written as 


dn 


Ni 


dt 


— ^v/sD^T (r<S'/i) 


-i (K+l f/\_ A 


SD 


ar 


K 


( 21 ) 


where /go is the fraction of newly born loops which die through the SD process, x ~ 0.5 is a 
numerical factor that characterizes the scaling configuration of the string network, Y ~ 100 
is a numerical factor that determines the life time of a loop due to gravitational radiation 
emission, and K ~ 0(1) is a numerical factor that determines the average length of the 
closed loops at their birth. It is generally argued G that ,/sd - h 

Using m, the above rate then gives the injection rate of massive Majorana neutrinos 
from cosmic string closed loops disappearing through the SD process per comoving volume 
as 

M l \ 

GeV/ ’ 


dY st 

al jVi 

dz / SD Z* \Vb-l 


1.57 x 10" 17 A7 / Mi 
-;- N N ‘ 


( 22 ) 


where z = Mi/T is the dimensionless variable with respect to which the evolution of the 
various quantities is studied. In the above we have used the following numerical values for 
the various constants: Y = 100, x = 0.5, g * = 100 and K = 1. 










2. Quick death 


Some fraction of closed loops may be born in configurations with waves of high har¬ 
monic number. Such string loops have been shown Q] to have a high probability of self- 
intersecting. In this case a loop of length L can break up into a debris of tiny loops of size 
(at which point they turn into the constituent massive particles) on a time-scale ~ L 
(hence “quick”). Since gravitational energy loss occurs over a time scale much larger than 
L, these loops do not radiate any significant amount of energy in gravitational radiation, 
and thus almost the entire original energy of these loops would eventually come out in the 
form of massive particles. 

Assuming again that each tiny loop of size ~ Vb-l yields a number Njy of heavy neutrinos, 
the injection rate of the massive Majorana neutrinos due to this “quick death” (QD) process 


is given by 


28] 


dn 


N i 


dt 


QD 


_ f AT 1 VB—L 


(23) 


where /qd is the fraction of loops that undergo QD. 

Note that while each tiny loop yields the same number Nn of the heavy neutrinos irre¬ 
spective of the SD or QD nature of the demise of the loop, the total injection rates of the 
heavy neutrinos in the two cases are different because of the different number of tiny loops 
that result from an initial big loop and the different time scales involved in the two cases. 

Using ED the above rate (T?0l) can be rewritten as 


dY st 

al Nx 

dz 


1.36 x 10” 36 


/qdA 


N 


QD 


VB-L \ S' Mi 

GeV / V GeV 


(24) 


While the value of /qd is not known, there are constraints [23, 41] on /qd from the 
measured flux of the ultrahigh energy cosmic rays above ~ 5 x 10 19 eV and more stringently 
from the cosmic gamma ray background (CGRB) in the 10 MeV - 100 GeV energy region 


measured by the EGRET experiment 


42]. The latter gives 


/qd?7i6 ~ 9.6 x 10 


(25) 


where rj i6 = (?1b-l/10 16 GeV). For GUT scale cosmic strings with ?/ 16 = 1, for exam¬ 
ple, the above constraint implies that /qd < 10 5 , so that most loops should be in non- 
selfintersecting configurations, consistent with our assumption of /sd ~ 1- Note, however, 
that /qd is not constrained by ( 1251 ) for cosmic strings formed at scales rj B _ L < 3.1 x 10 13 GeV. 










IV. LEPTOGENESIS IN PRESENCE OF COSMIC STRINGS: TOWARDS A 
COMPLETE ANALYSIS 

A. Analytical estimate 

In our previous work j^] we made simple analytical estimates of the maximum possible 
contributions of cosmic string loops to leptogenesis. Neglecting the contribution of the 
thermal Nr s, allowing the maximal value of the CP-violation parameter given by equation 
m and demanding that the resulting value of the baryon-to-photon ratio not exceed the 
observed value given by equation m we derived upper bounds on the mass M\ for the SD 
and QD processes of decay of cosmic string loops. This is shown in Figure [I] for N^ = 10 
for illustration. 



hB-L (GeV) 


FIG. 1: Constraints on hi = Mi/t/b-l from consideration of possible maximal contributions of 
B — L cosmic string loops to P-asymmetry. Models lying on the thick (thin) solid line can in 
principle produce the observed BAU entirely due to the SD (QD) process of B — L cosmic string 
loops. Models above the thick solid line are ruled out from consideration of overproduction of the 
B asymmetry. The h\ = 1 line is also shown for comparison. 







Figure UJ allows us to identify regions in the t)b-l — Mi plane determining models for 
which cosmic string processes discussed above can play a significant role in the leptogenesis 
mechanism. Clearly, from Figure El cosmic string are relevant for leptogenesis only for 
r) b-l A 1-7 x 10 11 GeV; lower values of t) B -l are relevant only if we allow h\ > 1. Note 
also that the above lower limit on t]b-l (for cosmic strings to be relevant for leptogenesis) is 
determined by the SD process; the QD process becomes relevant only at much higher values 
of TjB—Li namely, rf^_ L > 1.2 x 10 14 GeV. 

In this context, it may be mentioned here that, following our previous work, a recent 
work 113] has found that in the case of degenerate neutrinos (as opposed to hierarchical 
neutrinos assumed in our work), the B — L cosmic strings become relevant for leptogenesis 
at much higher values of t]b-l compared to those in the case of hierarchical neutrinos, e.g., 
t)b—l A 3.3 x 10 15 GeV for the SD process and i]b-l > 10 16 GeV for the QD process. As 
already mentioned, we assume hierarchical neutrino masses in the present paper. 

We now proceed towards a complete analysis of leptogenesis in presence of cosmic strings 
by first setting up and then solving the relevant Boltzmann equations that include the non- 
thermal Nrs produced from the decaying cosmic string loops as well as the already existing 
thermal Nr 1 s and also include all the relevant wash-out effects. 


B. Boltzmann Equations 


The total rate of change of the abundance of Ni s can be written as 


dY Nl 

dz 



+ 

D,S 



injection 


(26) 


The first term on the right hand side of equation © is given by the usual Boltzmann 
equation 


(^r) D =- (D + S) ( y "'- y n)’ (27) 

where D and S constitute the decay and A L = 1 lepton number violating scatterings, 
respectively, which reduce the number density of A 1; and is the abundance of in 
thermal equilibrium. 

The second term on the right hand side of equation © gives the rate of injection of 






Nfs due to disappearance of B — L cosmic string loops: 


dY, 


N i 


dz 


(dY&\ 

V dz ), 


+ 


dYft 

dz 


(28) 


QD 


and 


injection 

where the two terms on the right hand side are given by the equations 
respectively. 

The two terms in equation (121 compete with each other. While the first term reduces 
the density of Ah’s, the second term tries to increase it through continuous injection of Nf s 
from the disappearance of cosmic string loops. The CP-violating decays of the Nf s produce 
a net B — L asymmetry. This can be calculated by solving the Boltzmann equation 

dY B L 


dz 


= —e 1 D (Y Nl - Y£) - WY b _ 


L- 


(29) 


The first term in equation (El involving the decay term D produces an asymmetry while a 
part of it gets erased by the wash out terms represented by W which includes the processes 
of inverse decay as well as the A L = 1 and A L = 2 lepton number violating scatterings. In a 
recent work J44| it has been claimed that thermal corrections to the above processes as well 
as processes involving gauge bosons are important for the final B — L asymmetry. However, 
this is under debate j45|, and we have not included these corrections in the present paper. 

In equations El and El, D = T d/Hz, where 

1 


r D = 


is the decay rate of Ni, with 


167 TV 2 


(m/ D m D ) n 

m i = 


mi Mr 


the effective neutrino mass parameter 


Q. 


Mi 


(30) 


(31) 


and H is the Hubble expansion parameter. Also, 
S = T s/Hz, where Ts is the rate of A L = 1 lepton number violating scatterings and 
W = Yw/Hz , where Tw is the rate of wash out effects involving the A L = 1 and A L = 2 
lepton number violating processes and inverse decay. The various Ts are related to the 
scattering densities j^| y’s as 

(32) 


r?(z) = ^ z 


n 


eq 

X 


The dependence of the scattering rates involved in A L = 1 lepton number violating processes 
on rhi and Mi are similar to that of the decay rate T d■ As the Universe expands these T’s 









compete with the Hubble expansion parameter. In a comoving volume we have (with same 
notations as in |l|) 


7 d 


JVi 
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(33) 


sH(Mi) J ’ \sH(Mi) I ’ \sH(Mi) 

On the other hand the 7 ’s for the AL = 2 lepton number violating processes depend on mj 
and Mi as 


7tvi 


r Nl , t 


oc k^rfi^Mi. 


(34) 


sH(Mi) J ’ l sM(Mi) 

In the above equations (E3D, (ED, ki, i = 1, 2, 3 are dimensionful constants determined from 
other parameters; see |9|] for details. 


C. Constraint on the effective neutrino mass (mi) 

In order to satisfy the out of equilibrium condition the decay rate of Ni has to be less 
than the Hubble expansion parameter. This imposes a constraint on the effective neutrino 
mass parameter rh 1 as follows. At an epoch T < Mi, 


-jr = — — K < 1 , 

H 771* 

where rri *, the cosmological neutrino mass parameter , is given by 

^-yl /2 

m* = — = 6-4 x lCT 4 eV, 

V2G f 


(35) 


(36) 


where g* is the effective number of relativistic degrees of freedom. We see that a net B — L 
asymmetry can be generated dynamically provided rhi < m* at an epoch T < Mi. This 
constraint on rhi can be realized in a model as follows. 

We assume a charge-neutral lepton symmetry for which we take the texture of the Dirac 

3 

/ „ „ \ 

^m e m T ■ (37) 

0 y/m e m T m 7 


mass of the neutrino of the form 


m D = 


V 


/ 


Using m and ED in equation ED we get the constraint to be 

rhi < 5.25 x 10" 6 eV, 

which is in consonance with the requirement that mi < m* at any epoch T < M\. 


( 38 ) 















D. Solving the Boltzmann equations 


At an epoch T 3> M\ the lepton number violating processes are sufficiently fast as to set 
the B — L asymmetry to zero. As the temperature falls and becomes comparable to Mj, 
a net B — L asymmetry is generated through the CP- violating decay of N\. The resulting 
asymmetry can be obtained by solving the Boltzmann equations. We solve equations © 
and © numerically with the following initial conditions 

= an d Yj?_L = o. (39) 

Using the first initial condition we solve equation (|2Hj) for Yj ^ 1 , and the corresponding B — L 
asymmetry Yb-l is obtained from equation © by using the second initial condition of 
equation ©. 

In the usual thermal scenario the B — L asymmetry depends on fh\ and M\. In presence 
of cosmic strings there is in addition an explicit dependence on t)b~l since the injection 
rate of NrS from the cosmic string loops explicitly depends on it. To illustrate the effect 
of the cosmic string generated non-thermal ATr’s on the evolution of the B — L asymmetry, 
we display in Figures [2] and d] the results obtained by numerically solving the Boltzmann 
equations described above for some specific values of the various relevant parameters. It 
can be seen from Figure [21(a), for example, that in the absence of cosmic strings the B — L 
asymmetry approaches the final value at T ~ O.lMi when all the wash out processes fall out 
of equilibrium. In contrast, in the presence of cosmic strings, for sufficiently large values of 
Mi, the B — L asymmetry continues to build up until the injection rate of Nr s from cosmic 
string loops becomes insignificant. For h\ = 1 this happens around T ~ M\f 600 (see Figure 
lib)), which is far lower than in the purely thermal case. As a result of this, in presence 
of cosmic strings, for a fixed value of the symmetry breaking scale t\b~l = 10 13 GeV, for 
example, the final B — L asymmetry is enhanced by three orders of magnitude for the effective 
neutrino mass rhi = 10 -4 eV (Fig. Hb)) and by two orders of magnitude for rhi = lCU 5 eV 
(Fig. Hb)). For rjs-L = 10 13 GeV, the effect of cosmic string essentially disappears for 
hi < 0.01 

The above results, in particular the dependence on the value of m 1? can be understood as 
follows: First let us consider the purely thermal leptogenesis case, i.e., in absence of cosmic 
strings. Below a certain temperature when the wash out processes fall out of equilibrium 
the asymmetry produced by the decay of Ah’s does not get wiped out, and the produced 
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FIG. 2: Evolution of the B — L asymmetry for rhi = ICG 4 eV, tjb-l = 10 13 GeV and different 
values of M\ (a) in absence of cosmic strings, and (b) in presence of cosmic strings. The CP 
violation parameter e\ has been given its maximal value. 

B — L asymmetry remains as the final asymmetry. Since the decay rate of N\ s depends 
linearly on fh\ (see equation ©), a larger value of fhi implies that the condition of decay in 
out-of-equilibrium situation is satisfied at a later time when the abundance of thermal Ah’s 
is smaller, thus yielding a smaller final value of Yb~l■ And as expected, this effect is larger 
for larger values of M\. Now, in the presence of cosmic strings, for values of i]b-l for which 
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FIG. 3: Evolution of the B — L asymmetry for mi = ICG 5 eV, tjb-l = 10 13 GeV and different 
values of M\ (a) in absence of cosmic strings, and (b) in presence of cosmic strings. The CP 
violation parameter ei has been given its maximal value. 

the contribution of the cosmic string generated Ah’s to Y B _ L always remains insignificant 
compared to that due to the thermal Ah’s, the dependence of the final value of Y B _ L on ifii 
is essentially the same as in the absence of strings as explained above. However, for those 
values of t]b-l for which the string contributions to Yb~l are significant compared to the 
thermal contribution, the dependence of the final value of Yb-l on fh\ is opposite to that in 




























the absence of cosmic strings, i.e., the string contribution increases with increasing values 
of 777.1- This is easy to understand: The string generated Nr s make dominant contribution 
to Yr_l, if at all, only at late times (T -C Mi) when the abundance of the thermal Nr s 
has fallen to insignificant levels. The decay of the string generated Nr’s at such late times 
automatically satisfies the out-of-equilibrium condition. In this situation a larger value of 
rhi simply implies a larger rate of decay of the Nr ’s leading to a quicker development of the 
string contribution to Yr~l ■ This, together with the fact that the injection rate of the Nr s 
from cosmic string loops is higher at earlier times, leads to a higher final value of Y B -l- 
Figures Efb) and 01(b) plotted for rhi = 10~ 4 eV and fh\ = lCU 5 eV simultaneously bear out 
these expectations. 

Note that in Figures El and 01 the CP asymmetry parameter has been taken to have its 
maximum value, i.e., sin 6 = 1 in equation m ■ Clearly, for this value of the produced B 
asymmetry for i/b-l = 10 13 GeV and Mi >, 10 10 GeV (i.e., h\ >, 10 -3 ) exceeds the observed 
value given by Y§^ L ~ O(10 -10 ) even for purely thermal lcptogenesis. Including the effect 
of cosmic strings only makes the situation worse. This implies that we must place an upper 
bound on the CP asymmetry phase 5 depending on the values of 7]b-l : Mi and rhi, hr order 
not to overproduce the B asymmetry. At the same time, to produce the observed BAU 
for t]b-l = 10 13 GeV and hi = 1, for example, we require sin 5 = sin<5 req d = 0(1CT 7 ) and 
sin<5 req d = O(10 -6 ) for rhi — 10 _4 eV and rhi — 10 _5 eV, respectively. These values of sin 6 
diminish the purely thermal contributions, Yj £_ L , to 0(1CT 13 ) and 0(1CT 12 ), respectively. 

We have calculated the value of Y b ~l for a range of values of 7j B -L , hi and rhi, both 
in the purely thermal case, Y$_ L , and with the effect of cosmic strings included, Y^^. 
The results for r} B -L = 10 12 GeV and 7) B ~l = 10 11 GeV are summarized in Tables |H and in 
respectively, where we also indicate the order of magnitude of sin (5 req d, the value of sin 5 
required to produce the observed B asymmetry, for each set of the chosen parameters. 

From these Tables as well as from Figures [2] and 01 we see that there exist ranges of values 
of the parameters t/b-l, hi, fri\ and 5 for which the cosmic string generated Nr’s make the 
dominant contribution to and indeed produce the observed BAU while the purely thermal 
lcptogenesis mechanism is not sufficient. 


hi 

rhi = 10 4 eV 

rhi = 10 5 eV 

1 B-L 

-f^th+st 

1 B-L 

log (sin (5 req d) 

xrth 

1 B-L 

■^th+st 

1 B-L 

log (sin d reqd ) 

1 

2.72 x 10 -7 

1.23 x 10 -5 

-5 

3.90 x 10~ 7 

1.59 x 10~ 6 

-4 

0.1 

3.10 x 10“ 8 

3.22 x 10“ 8 

-2 

3.97 x 10^ 8 

3.98 x 10" 8 

-2 

0.01 

3.25 x 10“ 9 

3.25 x 10“ 9 

-1 

4.04 x 10^ 9 

4.04 x 10“ 9 

-1 

0.001 

3.39 x 10 _1 ° 

3.39 x 10~ 10 

0 

4.11 x 10~ 10 

4.11 x 10“ 10 

0 


TABLE I: The B — L asymmetry in the purely thermal leptogenesis case (Y^_ L ) and in presence 
of B — L cosmic strings for ijb-l = 10 12 GeV and different values of the Yukawa coupling 

h± and the effective neutrino mass rh\. The order of magnitude of the value of sin <5 required to 
produce the observed BAU, sinUeqd, for each set of values of the parameters fh\ and hi, is also 
given. 


hi 

rhi = 10 4 eV 

rhi = 10 5 eV 

Vth 

1 B-L 

■u^th+st 

1 B-L 

log(sin <5 reqd ) 

1 B-L 

■^th+st 

1 B-L 

log(sinh reqd ) 

1 

3.10 x 10~ 8 

1.51 x 10“ 7 

-3 

3.97 x 10 -8 

5.17 x 10~ 8 

-2 

0.1 

3.25 x 10“ 9 

3.27 x 10“ 9 

-1 

4.04 x 10“ 8 

4.04 x 10~ 9 

-1 

0.01 

3.39 x 10 _1 ° 

3.39 x 10 _1 ° 

0 

4.11 x 10“ 10 

4.11 x IQ’ 10 

0 


TABLE II: Same as Table H) but for tjb-l = 10 11 GeV. 

V. CONCLUSION 

We have studied the effect of B — L cosmic strings arising from the breaking of a U(1)b-l 
gauge symmetry, on the baryon asymmetry of the Universe. The disappearance of closed 
loops of B — L cosmic strings can produce heavy right handed neutrinos, Nr s, whose 
CP-asymmetric decay in out-of-thermal equilibrium condition can give rise to a net lepton 
(L) asymmetry which is then converted, due to sphaleron transitions, to a Baryon (P) 
asymmetry. We have solved the relevant Boltzmann equations that include the effects of 
both thermal and string generated non-thermal N^s. By exploring the parameter region 

























spanned by the effective light neutrino mass parameter mi, the mass M\ of the lightest of 
the heavy right-handed neutrinos (or equivalently the Yukawa coupling h\) and the scale of 
the B — L symmetry breaking, tjb-lj we found that there exist ranges of values of these 
parameters, in particular with t)b-l > 10 11 GeV and h\ >, 0.01, for which the cosmic string 
generated non-thermal Nr s can give the dominant contribution to, and indeed produce, 
the observed Baryon Asymmetry of the Universe when the purely thermal leptogenesis 
mechanism is not sufficient. We have also discussed how, depending on the values of t)b—Li 
rhi and hi, our results lead to upper bounds on the CP violating phase 5 that determines 
the relevant CP asymmetry in the decay of the heavy right handed neutrino responsible for 
generating the L-asymmetry. 
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